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ABSTRACT 

We analyse the effect that peculiar velocities have on the cosmological inferences we make using 
luminosity distance indicators, such as type la supernovae. In particular we study the corrections 
required to account for (a) our own motion, (b) correlations in galaxy motions, and (c) a possible 
local under- or over-density. For all of these effects we present a case study showing the impact on 
the cosmology derived by the Sloan Digital Sky Survey-II Supernova Survey (SDSS-II SN Survey). 

Correcting supernova redshifts for the cosmic microwave background (CMB) dipole slightly over- 
corrects nearby supernovae that share some of our local motion. We show that while neglecting the 
CMB dipole would cause a shift in the derived equation of state of Aw ~ 0.04 (at fixed f2 m ) the 
additional local-motion correction is currently negligible (Aw < 0.01). 

We then demonstrate a covariance-matrix approach to statistically account for correlated peculiar 
velocities. This down-weights nearby supernovae and effectively acts as a graduated version of the 
usual sharp low-redshift cut. Neglecting coherent velocities in the current sample causes a systematic 
shift of Aw ~ 0.02. This will therefore have to be considered carefully when future surveys aim 
for percent-level accuracy and we recommend our statistical approach to down-weighting peculiar 
velocities as a more robust option than a sharp low-redshift cut. 
Subject headings: cosmology: observations — supernovae : general 
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Concordance cosmology (ACDM) is a successful 
model of ou r universe, fitting observations of type la 
supernovae (IRiess et al.l 119981; iPerlmutter et alT 119991: 



Astier et all 120061: iWood-Vasev et all 120071: IRiess et al l 
i 2004 120071: iKowalski et all 120081: jHicken et all 120091: 
iKessler et all 120091 : iFreedman et all 2009), the cosmic 
microwave background fPage et alJl2003t iTeemark et al 
ll20f) 



[20061: ISnergel et all [2007; iKomatsu et all 120091 



2011 



baryon acous tic osc illatio n s (lEisenstein et^T 120051 : 



iPercival et all [20071 [20Tot iBlake et all l2011blh and 
growth of structure (|Blake et all 1201 lal) . amongst oth- 
ers. However, it relies on the existence of dark compo- 
nents of the universe - dark energy and dark matter - 
whose nature remains mysterious. This has given rise 
to questions about the validity of our theory of gravity 
itself. Therefore enormous observational effort is contin- 
uing to better characterise the dark sector by measuring 
the expansion history of the universe and the growth of 
structure within it. 

Type la supernovae remain a lynch-pin in this effort, 
and more surveys are underway, or planned, to gather 
ever more high quality data to try to reduce the uncer- 
tainties on our cosmological parameters down below the 
1% level and search for possible variations in the equa- 
tion of state of dark energy. To achieve this accuracy 
we will have to address small systematic effects that had 
previously been negligible. In this paper we consider sys- 
tematic errors that could arise from neglecting the pecu- 
liar velocities and gravitational redshifts induced by large 
scale structure. 

The customary diagnostic in supernova cosmology is 
the Hubble diagram, a measurement of luminosity as 
a function of redshift. When using type la supernovae 
to measure this magnitude-redshift relation, the redshift 
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used should be entirely due to the expansion of the uni- 
verse. In practice this is never the case, as large-scale 
structure in the universe induces peculiar motions so 
that the measured redshift contains some contribution 
from peculiar velocities. To date, the majority of the ef- 
fort in calibrating type la supernova measurements has 
been increasing the accuracy and precision with which 
we can determine their luminosity, and thus their use as 
a standard candle. In comparison the uncertainty on the 
redshift of the supernovae has usually been considered 
negligible. It is this more neglected uncertainty we turn 
our attention to in this study. 

We distinguish between the cosmological redshift, z 
due entirely due to the expansion of the universe; the 
'peculiar' redshift, z pcc , due entirely to peculiar veloci- 
ties; and the gravitational redshift, z gra v, due to density 
fluctuations. These three redshifts combine to give the 
observed redshift, z, according to, 

(1 + z) = (1 + z)(l + z pcc )(l + z grav ). (1) 

Systematic peculiar velocity effects can manifest them- 
selves in a number of ways. In analogy to the cosmic mi- 
crowave background (CMB), we can consider monopole, 
dipole, and higher-multipole effects. 

• Monopole effects mimic expansion by imprinting 
a non-cosmological redshift isotropically on all 
sources. They can arise if we are in the center of a 
local over- (or under-) density due to both peculiar 
velocities of infall (or outfall) and gravitational red- 
shifts. Due to their isotropic nature they are very 
difficult to distinguish from cosmological redshifts 
and could lead to significant systematic errors in 
cosmological parameter estimation. 

• Dipole effects occur due to our own peculiar mo- 
tion, or due to being off-center in a local over- 
or under-density. Our motion with respect to the 
overall expansion is well measured by the dipole in 
the CMB, and can be easily corrected-for. How- 
ever, the correction needs to be modified when cal- 
culating our velocity relative to local galaxies that 
share some of our motion with respect to the CMB. 

• Higher-order multipole effects occur when distant 
sources move coherently as they share gravitational 
attraction to common large-scale structure. These 
are analogous to the fluctuations seen in the CMB 
after the monopole and dipole terms have been re- 
moved. 

All these peculiar velocities are manifestations of coher- 
ent flows in the universe. We consider the lowest order 
coherent flows (monopole and dipole) separately because 
they have the largest systematic effects. 

When we compare the precision of redshift (z) mea- 
surements with the precision of apparent magnitude (m) 
measurements we are primarily concerned with their ef- 
fect on the Hubble diagram (rather than the relative pre- 
cision as measured by Az/z and Am/m). The effect on 
the Hubble diagram can be quantified by considering the 
slope of the magnitude-redshift relation, dm/dz. Indeed, 
the standard method for including redshift uncertainties 
(a z ) in cosmological analyses is to convert them to mag- 
nitude uncertainties (a m ) using the dm/dz derived from 



a fiducial cosmological model (see Appendix [S] for an 
outline of the procedure) . The conversion between a red- 
shift uncertainty and a magnitude uncertainty is shown 
in Figure [TJ for a few different redshifts. At higher red- 
shifts the magnitude-redshift relation is flatter {dm/dz is 
smaller) which means large redshift uncertainties gener- 
ate only small magnitude uncertainties. At low redshifts 
the converse is true, and small redshift uncertainties give 
large magnitude uncertainties. 

The measurement of redshift will remain far more accu- 
rate than the measurement of the supernova magnitude 
into the foreseeable future. However, the accuracy of 
those measurements can be misleading since systematic 
effects on redshift due to peculiar motions can be much 
larger than the measurement error (see Fig. [I}. 

The dominant source of intrinsic redshift dispersion 
(that is always included in cosmological analyses) is the 
effect of random peculiar velocities. These are usually 
taken to be abou10 <jP cc = 300km s" 1 , which accord- 
ing to w pcc = cz pcc corresponds to an error in redshift 
of cr poc = 0.001. This redshift uncertainty gives a non- 
negligible magnitude uncertainty of <jf£ c = 0.2 for ob- 
jects at redshift z = 0.01, which reduces to a^ c — 0.02 
for objects at z = 0.1. These values should be compared 
with the intrinsic diversity in supernova magnitudes of 
a rn ~ 0-1 an d the observational magnitude uncertainty 
(including the uncertainty in fitting the SN light curves) 
of tr™ cas < 0.1 for the most distant supernovae included 
in most samples (closer ones have smaller observational 
uncertainty) . 

The observational uncertainty in redshift is typically 
fi| pcc ~ 0.0005 for SN host-galaxy-based redshifts mea- 
sured by SDSS, and cr| poc w 0.005 for redshifts based on 
the supernova spectra alone (jZheng et al.l 120081 ). So the 
observational uncertainty is of similar magnitude to the 
intrinsic scatter due to peculiar motions. 

Treating these motions as random scatter is not en- 
tirely valid, since galaxies (and the supernovae in them) 
preferentially fall into overdense regions, so objects in 
the same region of sky will tend to have correlated pe- 
culiar velocities. Detailed measurements of this effect 
may provide another technique for measuring the matter 
distribution of the universe and thus derive cosmological 
parameters using diagno stics such as the peculiar veloc- 
ity power spectrum ( e.g. iBonvin et alJl2006aJ: iNeill et all 
20071 iGordon et all 120071 120081: lAbate k Lahavl 120081: 
Hannestad et all 120081 ) . iLampeitl et all (|2010l App. B) 
showed that the signal in the SDSS SN data set is 
too small to measure cosmological parameters this way 
(mostly because our sample is too distant). Here we 
therefore concentrate only on the deleterious impact pe- 
culiar velocities have on the cosmological results derived 
from the SN magnitude-redshift relation. 

1.1. Data 

In this paper we use the Sloan Digit ial Sky Suryey-II 
Supernova Survey (SDSS-II SN Survey lYoidTeFall l2000l: 

18 For the SDSS sample IKessler et all J2u09t) use o-g ec = 
300km s — 1 for the random peculiar motions added in quadra- 
ture with 0"5 ec = 200km s _1 for the internal velocities, giv- 
ing a total crP cc =360km s -1 , corresponding to cr^ ac of 0.0012. 
For the ESSENCE sample IWood-Vasev et al.l 12007T ) uses cr% cc = 
400km s" 1 . 
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TABLE 1 

Sub-samples of type Ia supernovae used by Kessler et al. 
2009. The top three rows give the survey name, redshift 
range, and number of supernovae in each data set. below 
the line, ticks indicate which data sets are included in 
each of the sub-samples. 
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Holtzman et~aT1 120081 : iFrieman et all l2008t iSako et al.l 
2008J) as a "case study. For the SDSS-II SN Survey 
repeat images were taken of an equatorial stripe, 2.5 
degrees wide and about 120 degrees long centered on 
RA 23.5hr (SDSS Stripe 82), which corresponds to a 
field center at Galactic coordinates (£,b) « (84°,— 57°). 
The direction of the CMB dipole in Galactic coordi- 
nates is toward (e,b) = (263°.85 ± 0°. 1,48°. 25 ± 0°.04) 
([Bennett et al.1 12003T ); so the antipode direction lies at 
(£,b) « (83°. 85, -48°.25). That means that the center 
of the SDSS field is almost aligned with the direction in 
which the CMB appears blue-shifted. 

Over three years about 500 spectroscopically confirmed 
Type Ia supernovae were discovered in the redshift range 
0.05 < z < 0.4. The first y ear's data includ i ng 10 3 su- 
pernovae w ere published i n iHoltzman et al.l (|2008fi and 
analysed in iKessler et al.l (|2009t hereafter referred to as 
K09 j3 who combined these with a re-analysis of existing 
data to make a coherent sample of 288 supernovae that 
were used to measure cosmological parameters. It is this 
data set that we use here, focussing on several of the 
subsets they define to demonstrate the effect of different 
redshift ranges (see Table [I}. 

Since the SDSS sample concentrates on relatively 
nearby supernovae (0.05 ^ z < 0.4) the peculiar velocity 
contribution is a more significant proportion of the to- 
tal redshift than in surv eys that focus on higher red- 
shifts such as ESSENC E (IWood-Vasev eFaTI f2007h . Su- 
pe rNova Legacy Sur vey (jAstier et al.ll2006[ ). and Higher- 
z (|Riess et al.l 120071) . Moreover, since the region of sky 
covered by the SDSS SN survey lies close to the direction 
of the CMB dipole, (Fig. [5]) the alignment conspires to 
maximise the magnitude of the effect. For all these rea- 
sons the SDSS SN sample is an interesting one in which 
to test peculiar velocity effects. 



1.2. Notation 



We use the metric 
ds 2 = -c 2 dt 2 



R(t) 2 [d X 2 + S 2 ( X )d® 2 



(2) 



where t is proper time, R(t) is the scale factor with 
dimensions of length, \ is the dimensionless comov- 
ing coordinate, Sk(x) — sm (x)i Xi sm h(x) for closed, 
flat, and open universes respectively, and d® encom- 
passes the angular terms. The present day scale factor, 
R Q = c/(-ff \/|f2 fc |), an d the dimensionless scale factor is 

19 http: //das . sdss . org/va/SNcosmology/sncosm09_f its . tar . gz 
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Fig. 1. — Plot of the conversion between a redshift uncertainty 
a z and a magnitude uncertainty a m , for a variety of redshifts be- 
tween z = 0.001 and z = 0.02 as labelled. Dotted lines depict the 
typical redshift dispersion from random peculiar velocities, and the 
observational uncertainties in magnitude and redshift (the obser- 
vational uncertainty for cr z assumes we have redshifts from host 
galaxy spectra, when only SN spectra are available the redshift 
error is an order of magnitude higher). Those random contribu- 
tions are all marked '(rand)', while the other dotted line depict- 
ing 'CMB (syst)' shows the maximum systematic shift in redshift 
caused by the CMB dipole. This figure shows why the low-redshift 
cutoff of z = 0.02 is appropriate, because above this redshift all 
redshift uncertainties are significantly smaller than the observa- 
tional magnitude uncertainty. The conversion has been done for 
the (n m ,f2/v) = (0.3,0.7) model, but there is negligible difference 
between this and most plausible homogeneous cosmological models 
over this redshift range (including the empty model and a model 
with f! m = 0.3 but no cosmological constant). 

defined as a = R/Rq. Hubble's parameter is H(z) = a/ a, 
where an overdot represents differentiation with respect 
to proper time. The dimensionless comoving distance as 
a function of redshift is related to cosmological parame- 
ters by 

X=^[^- (3) 



Rq 



H( z y 



Frequently we will need the comoving distance with units 
of length, for which we will use the shorthand \ = RoX- 
We begin in Section [2] by providing the theory needed 
to make dipole and correlated velocity corrections to dis- 
tance measurements. Although focusing on type Ia su- 
pernova measurements this analysis is general to any lu- 
minosity distance measure. In Section[3]we then perform 
a case study of these effects on the cosmological infer- 
ences made using the supernova data from the SDSS II 
supernova survey (Kes sler et al.ll2009h . We conclude in 
Section [5] 

2. CORRECTING FOR PECULIAR VELOCITIES 

2.1. Dipole correction 

When using supernova redshifts to make cosmologi- 
cal inferences we need to remove the imprint of our 
own peculiar motion so that the redshift of the super- 
nova is entirely due to the expansion of the universe. 
To first order this is straightforward, since we know our 
own velocity to high precision from measurements of the 
CMB dipole. Correcting for the CMB dipole is stan- 
dard practise in a l l supernova cosmology analyses (e.g. 
Astier et al.l 120061: IWood-Vasev et all 120071: iRiess et al.l 



20071 iKowalski et al. 1 120081: IKessleT et al. How- 
ever, no correction is typically made for our motion rel- 
ative to nearby galaxies, for which the CMB dipole is a 
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poor approximation. Here we investigate the impact of 
both effects. 

We are moving at v© ec = 371 kms -1 with respect to 
the CMB (jKogut et alJll99llBennett et all [20031) . This 
is small compared to the Hubble flow for all but the near- 
est objects, dropping to less than ~ 1% beyond a red- 
shift of 0.1. Our motion thus contributes a maximum 
redshift change of a z = 0.00124 to sources that are di- 
rectly aligned with the dipole. This is much less than 
the equivalent uncertainty in our magnitude measure- 
ment (see Fig. [T]); it is only its coherence among SNe 
that can make it significant. 

The general relationship between redshift and 'peculiar 
velocity is, 



pec 



l + Z 



which simplifies to z pc 



j l + vP cc /c 
1 - vP cc /c 



(4) 



pec 



I c in the non-relativistic 
limit The redshift correction required to account for 
our velocity with respect to the CMB, v©, is 



(-n)/c 



(5) 



where n is the unit vector from the sun to the supernova. 
(The negative sign ensures that if we are moving in the 
direction of the supernova the resulting correction is a 
blueshift.) 

The observed heliocentric redshift, z, is then related to 
the cosmological redshift, z, byJHH 



(l + z) = (l + z)(l + z ec ). 



(6) 



Note that the NED velo city calculator 

(|NASA/IPAC Extragalactic Database! [2001 ) uses 
the approximation, 

z^z + z^. (7) 

This gives a fractional erroi[3 of precisely z Q 0C , which is 
negligible for most circumstances. 

The dipole not only shifts the redshift but also changes 
the apparent magnitude of the source due to the Doppler 
shift of the photon energy and relativistic beaming. The 
CMB dipole therefore also has an effect on the luminosity 
distance calcu l ated from the magn i tude o f a supernova 
dSasakil 119871 iPvne fe Birkinshawl 119961: iBonvin et al.l 
I2006al iCoorav fe Caldwelll 120061 : iHui fe Greend 12006) . 
This arises because the luminosity distance is related to 
the comoving distance, x> by (recalling that overbars re- 
fer to observations made from the CMB rest frame), 



d L (z) = (l + z)R Q S k ( X )- 



(8) 



However, what we actually observe is (recalling that z 
is the observed redshift and considering for the moment 

20 Note that the special relativistic velocity-redshift relation in 
Eq. |4]is only appropriate for peculiar velocities. It is not appropri- 
ate for recession velocities (the velocity that appears in Bubble's 
Law) , for which special relativistic cor rections should never he ap- 
plied IjDavis fc Lineweaver 2004, 2003). 

21 This assumes the observer has already corrected both for 
the motion of the Earth around the Sun, which contributes up 
to 30km s -1 depending on the time of observation and for atmo- 
spheric refraction, which contributes up to 90km s -1 (the index of 
refraction of air is 1.0003, so Az = 0.0003 and cAz = QOkms" 1 ). 
Usually this is done as a standard step in wavelength calibration. 

22 Rearranging Eq. [6] and Eq. [7] gives - — _ ~ z = 2q SC . 
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Fig. 2. — Map of the distribution of type la supernovae in the 
Nearby, SDSS, ESSENCE, SNLS, and HST samples in Galactic 
coordinates with the local (z ~ 0.1) dipole indicated by diamonds 
and the CMB dipole indicated by stars. The center of the SDSS-II 
SN stripe lies close to the direction of the CMB dipole. which makes 
it important to carefully correct the SDSS sample for the effects of 
the dipole. The local dipole is measured with respect to the CMB 
and is in approximately the opposite direction since galaxies in our 
local neighbourhood tend to share some of our peculiar velocity 
with respect to the CMB. 



only our own motion), 



d L (z) = (l + z)R Q S k (x), 

= (l + z)(l + z p Q cc )R S k (x) 
= J L (z)(l + z cc ). 



(9) 
(10) 

(11) 



So both the redshift and the luminosity distance need to 
be corrected for the effect of the dipole 

Alternatively, one can choose to correct for both in 
one fell swoop by correcting the observed luminosity dis- 
tance at redshift z to the luminosity distance that would 
have been observed at redshift z in the absence of pe- 
culiar velocities. This is the approach taken by the 
commonly used program simple_cosfitter (jConlev et al.l 
l2TJ0ll l^ lHi"u fc Greend (pOM) give the formula for d L (z), 
which can be used to correct only the g?l values without 
co rrecting z. Consideri ng only our own motion, Eq. 15 
of IHui fc Greenel (2006) can be rearranged to give, 



d\,{z) = d h (z) 



1 



a' e R a T k (x) 



vo-n 



(12) 



where a' e = da e /dr represents the derivative of the scale 
factor with respect to conformal time, evaluated at the 
time of emission^] and here we have kept the curva- 
ture dependence explicit, with T k (x) = tan(x),x, and 
tanh(x) in closed, flat, and open universes respectively. 

23 You may be concerned that in going from Eq. [9] to Eq. 1101 
we've neglected the factor of z in the calculation of \ = 
(c/Ro) Jq dz/H(z). However, this cosmological redshifting is in- 
dependent of the motion of the emitter or observer, and therefore 
does not need correcting for peculiar velocities. As long as we cor- 
rect the redshift of the supernova to the CMB frame our theoretical 
model comparison wil l be correct. 

24 http:/ /qold. astro. utoronto.ca/conley/simple_cosfitter/html/ 

25 We give conformal time dimensions of time, so dr = dt/a 
and the conformal time derivative is related to the proper time 
derivative (denoted by an overdot) according to 



da dt 
dt dr 



(13) 
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Fig. 3. — The velocity two-point correlation, £12, in units of 
(kms -1 ) 2 . The upper panel shows the correlation as a function 
of Z2 (the redshift of the higher redshift supernova) compared to 
a low-redshift supernova at a fixed z\ = 0.01, for three different 
angular separations (0, 30 and 100 degrees). As the redshift sepa- 
ration increases the correlation diminishes. The lower panel shows 
the correlation as a function of angular separation, in the case 
where the two supernovae are at the same redshift (2:2 = z\ =0.01, 
0.02, 0.05, and 0.1). For a fixed angular separation the correlation 
is most dramatic at low redshifts because this corresponds to a 
smaller physical distance than the same angular s epar ation at high 
redshifts. £12 is given by the last two lines of Eq. 1291 



When there are two peculiar velocities to correct, such 
as when accounting for the supernova's motioro with 
respect to the CMB (-Zg^ c ) m addition to our own motion, 
one uses 



(1 + z) = (1 + z)(1 + 4 cc )(1 + 4n)- 



(15) 

This equation is valid even for relativistic velocities, but 
in the literature it is mo re common to encoun ter approx- 
imate formulae such as (Hui & Greene 20061). 

(1 + z) = (1 + z)(l - v Q • n/c + v SN • n/c), (16) 
= (l + z)(l + Z l° c + z^ c ), (17) 

which are perfectly appropriate for the low velocities we 
encounter in almost all practical situations. 

When including the source motion the correction to 
the luminosity distance becomes, 



d L (z) = d L (z)(l + zP° c )(l + zP° c ) 2 



(18) 



Note that two factors of (1 + z<^) enter the luminosity 
distance correction. One is due to the Doppler shifting 
of the photons, the other is due to relativistic beaming. 

2.2. Coherent flows 

The next peculiar velocity effect we consider is the im- 
pact of coherent bulk motions. Large scale structure in- 
troduces correlated peculiar velocities as neighbouring 



26 The additional redshift due to the supernova's motion is 

7c = vg° c ■ n/c, (14) 
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where again n is the unit vector from the sun to the supernova and 
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Fi g. 4 . — The magnitude covariance due to peculiar motion 
(Eq. [29]). As for Fig. |3] the upper panel shows, for three different 
angular separations, the covariance between a z\ = 0.01 source 
and a source at a higher redshift, Z2- The lower panel shows 
for fixed redshifts (z 2 = z x =0.01, 0.02, 0.05, and 0.1) how the 
covariance drops as the angular separation increases. To get a 
feel for how much additional uncertainty the correlated compo- 
nent adds to the uncorrelated dispersion consider that for these 
four redshifts the canonical 300km s -1 dispersion corresponds to 
ct 2 ^ ~ (25,10,7,5) X 10~ 3 respectively (see Appendix [A} . So the 
strongest correlated uncertainties are an order of magnitude lower 
than the random dispersion. However, the random dispersion can 
be beaten down by raising the number of SNe, while the correlated 
covariance cannot. Thus, as the total number of SNe increases the 
correlated noise becomes comparatively more and more important. 



galaxies, and the supernovae they host, fall towards the 
same overdensities. Ignoring these correlations under- 
estimates the uncertainty in our cosmological inferences. 
The effect is particularly important at low redshift where 
supernovae will tend to be physically closer to each other 
(as a function of angular separation on the sky) . 

These coherent flows are the higher-order manifesta- 
tions of the gravitational influence of large scale struc- 
ture, beyond the local dipole we discussed in Sect. 12.11 
This exercise is analogous to subtracting the dipole from 
the CMB and considering the residual fluctuations. 

The reason correlated velocities have a deleterious im- 
pact on standard supernova cosmology arises because two 
correlated supernovae can not be statistically averaged in 
such a way so as to reduce the error in proportion to the 
square-root of the number of supernovae, as is usually 
assumed. Correlations mean the data are not randomly 
distributed about the central value, so the uncertainties 
do not average out exactly and there remains a resid- 
ual error even in the limit of an infinite number of data 
points. 

Here we calculate the expected statistical correlations 
in the peculiar velocities of galaxies as a function of their 
distance from each other based on a ACDM model. We 
convert this into an observationally useful measure by 
converting it to the expected covariance in magnitudes as 
a function of angular separation and redshift separation. 

Recall that supernova cosmology aims to fit the obser- 
vations of apparent magnitude, m, and redshift, z of a 
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supernova to the following relation, 

m(z) = 5log 10 d,L(z) + M, 



(19) 



where dh is the luminosity distance (in units of lOpc) and 
M is a constant incorporating the absolute magnitude 
of the supernova and Hubblc's constant. The luminosity 
distance is a function of the cosmological parameters we 
want to fit. 

The likelihood of a particular model, in a gaussian dis- 
tribution, is proportional to e~ x ' . Let rhi represent 
the ith measurement and rrii the corresponding model 
prediction. When all data points are independent, x 2 is 
given bvPl 

(rhi - rrii) 2 



X 



(20) 



Here Oi = a mi is the magnitude uncertainty on rrii . How- 
ever, when data points are correlated the more general 
form of x 2 is given by, 



X 



= ^2(rhi- miJCt/irhj -rrij), 

i,3 



where the covariance matrix, 



C > 



(SrhiSrhj) , 



(21) 



(22) 



quantifies how likely two supernovae are to have the 
same offset from the correct model. The factor Srhi = 
rhi ~ (fhi) designates how far the ith data point deviates 
from the mean of the observational data. Correlations in 
these deviations arise for several reasons in addition to 
the peculiar velocity affect we consider here, for example 
iBlondin et all (|201lD studies the covariance between dis- 
tance prediction errors due to magnitude error, peculiar 
velocities, and intrinsic covariance. 

The covariance matrix can be divided into the random 
component Ui and the correlated component, which we 
consider to be only due to peculiar velocities, Cjf, 



c 



vcl 



(23) 

As discussed in lHui fc Green e (2006), there are a wider 
variety of large scale structure induced fluctuations than 
are accounted for in Eq. 1231 For instance, lensing intro- 
duces correlated noise in addition to Poissonian fluctu- 
ations. There are also fluctuations due to gravitational 
redshift and the integrated Sachs- Wolfe effect. It can be 
shown that all these effects can be neglected for surveys 
of current practical interest (|Hui fc Greend [20061 ) . 

The random uncertainties, which contribute to the di- 
agonal part of Cij, include the intrinsic diversity in the 
supernova population, cr™ tr , the scatter due to measure- 
ment uncertainty, erf 103,8 , and the contribution from ran- 
dom peculiar velocities, aj el . Although random peculiar 
velocities primarily add dispersion in the apparent red- 
shift of the sources (the effect on the luminosity is smaller 
and is usually neglected), this is usually converted to a 
dispersion in magnitude and added in quadrature to the 
other magnitude uncertainties (see Appendix [A] for more 
detail), 

(24) 



inti" 1 



nicas2 



27 In this paragraph \ 2 represents the statistic, not the comoving 
coordinate squared. 



Since random peculiar velocities are taken into account 
by this diagonal term we set all diagonal terms C™ 1 = 0. 

Alternatively we could remove aj cl and reinstate them 
as the diagonal elements of Cy . 
The velocity correlation function is defined to be, 



Si 



<(v; -Xi)(vj •Xj)) 



(25) 



where Xi and xj represent the unit vectors pointing to- 
wards SNe i and j respectively, and Vi and vj represent 
the velocity vectors of each supernova's motion. 

The peculiar-motion-induced magnitude covariance is 
related to the velocity correlation function £ 2 vel by, 



C 



vcl 



clnlO 



1 _ ^iiL 
a 't Xi. 



1-^ 

a'j Xj 



£-vcl 



(26) 



shifts 



where c is the speed of light, x = RoX is the radial co- 
moving distance, a = R/ Rq is the normalised scale fac- 
tor, and the prime denotes the conformal time derivative. 
All quantities with a subscript i or j are to be evaluated 
at the redshift of the SN in question. For a non-flat uni- 
verse x -> RoTk(x)- 
A numerical code to compute both and 
for a pair of points at arbitrary red- 
and angular separation in the standard 
cosmological model of ACDM is available at 
http : //www. astro . Columbia. edu/^lhui/PairV. 
We illustrate the results of that code in Figs. |3] and |U 
and in what follows we explain the theory behind those 
calculations. 

To calculate the expected velocity correlation function 
given a theoretical model we need information about how 
structure grows. To first order this is given by the linear 
growth factor D(z) = 6(z)/S(0), where the overdensity 
S = (p — (p))/ (p) ■ As input we also use the mass power 
spectrum of density fluctuations observed at the present 
time -P(fc) z= 0j where k is the comoving wavenumber (in- 
verse distance). 

Using these we can estimate the distribution of pecu- 
liar velocities expected in a particular theoretical model. 
Concentrating for the moment only on the dispersion 
(the diagonal terms in the velocity correlation func- 
tio n), one finds the dis persion in peculiar velocities to 
be (IHui & Greendl200ll . 

< cl2 = Qf = D'{ Zl f I ^P(k) z= o, (27) 
which results in a dispersion in apparent magnitude of, 
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(Zj C 






1 - 




clnlO 




a 'i Xi. 



(28) 



In principle this dispersion is sensitive to non-linear fluc- 
tuations, but the velocity power spectrum weights larger- 
scale modes more than the density power spectrum does 
and we find that when using the linear mass power spec- 
trum for a ACDM model the resulting value for a Vi agrees 
with the canonical value of 300 km/s to better than 10% 
for all redshifts of interest. The off-diagonal components 
of Cij should be at least as well fit by linear theory since 
they are less sensitive to small-scale structure than a Vi . 
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Fig. 5. — Peculiar velocity correction required for the SDSS sam- 
ple, plotted as a function of right ascension (J2000). The SDSS 
SN sample spans an equatorial strip, and thus the peculiar veloc- 
ity correction is systematic with right ascension. The CMB dipole 
correction is shown in black, while the local dipole correction, which 
is more relevant for low-redshift sources, is shown as shaded points, 
with different shades representing different redshift ranges. 



The off-diagonal part of C ih given by C™ 1 in Eq.[H 
accounts for the effects of correlated peculiar flows. 
Expressing th i s qu antity in an observer-centric form 
iHui fc Greene! (|2006[ ) show that for a flat universe, 



: In 10 



m 

oo 



1 



dk 
2^ 



1 



(29) 



P(k) z=0 x 



where Vi is the Legendre polynomial, jg is the spherical 
Bessel function, and j' e is its derivative with respect to 
its argument. It is useful to note that j'Ax) — je-i — {£ + 
l)ji/x. This observer- centric form can be d erived from 
Eq. 22, D7, and D10 of IHui fc Green! (|2006l ). by setting 
the survey geometry to be two delta functions localized 
at the two SNe of interest. 

An alter native separa t ion-centric form fo r the same 
quantity is (|Gorskilll988HGordon et aLl l2007). 



i 2 r 



clnlO 
H(r) = DM 



(30) 



f)n(r) + [x i -x j -(x i .f)(x j .f)]E(r)] 

2jx(kry 



2ir z 



S(r) 



dk 

2^2 



jo(kr) 

ji(kr) 
kr 



kr 



where the comoving separation between the two SNe is 
given by r and f is the unit vector pointing along the 
separation. 

That Eq. (f2"9"|) and ([31?)) are equivalent is shown in Ap- 
pendix B. The separation-centric form is useful for fast 
computation, while the observer-centric form is more di- 
rectly linked to the observed velocity angular power spec- 
trum. Note that these two equations are only strictly 



Fig. 6. — Similar to Figure [5] but with the correction expressed 
as a percentage and plotted against redshift. The peculiar velocity 
correction taking into account the local dipole (shaded filled cir- 
cles) is compared to the CMB dipole corection (black filled circles). 
The difference between the two is shown as shaded diamonds. The 
dashed line shows the theoretical maximum correction that would 
be applied to an object directly aligned with the CMB dipole: a 
constant peculiar velocity correction of 371km s" 1 that decreases 
with redshift only because it is a decreasing fraction of the total 
redshift. The difference between the CMB and local dipole cor- 
rection is only significant for the closest supernovae in the sample, 
with a ~4% correction in redshift when z < 0.02 but a correction 
of less than 1% for supernovae with z iS 0.05. 



valid for a flat universe, since the derivation in Appendix 
B uses a plane-wave expansion that needs modification 
if the universe is not flat. 

We are interested in the implications of deviating from 
the common practise of assuming all velocities are un- 
correlated. In Section |3~21 we therefore test the impact of 
including C™ 1 in the covariance matrix of the uncertain- 
ties for the supernova sample used by K09. 



3. IMPACT ON COSMOLOGICAL PARAMETERS 

3.1. Impact of the local dipole 

It is inappropriate to apply the full CMB dipole cor- 
rection to nearby supernovae because galaxies in our 
local universe share some of our locally induced pe- 
culiar motion. For example, below redshifts of ^0.02 
we are strongly influenced by th e Great Attractor an d 
the Perseus-Pisces S upercluster (jErdogdu et al.l [20061) . 
iBonvin et al.l (|2006b[ ) measured our dipole r elative to the 
nearb y sample of 44 supernovae used by (jAstier et al.l 
I2006T ) and found it to be consistent with the CMB 
dipole, although with large uncertainties (about ±30° 
directio nal uncertainty, and 200 km s _ 1 magnitude uncer- 
tainty). IHaugbolle et al.l (|2007l ) were able to more pre- 
cisely measure the velocity flow of the local uni verse using 
the 13 3 low-redshift type la supernovae from Uha et aLl 
(|2007l) . At a redshift of - 0.02 (60/i- x Mpc) they find 
a dipole amplitude of 239lgg kms _1 in the direction 
£ « 281° ± 23°, b w 14° ± 16° (measured relative to 
the CMB rest frame). In Figure [5] this dipole is marked 
by diamonds. The magnitude of this dipole decreases 
with redshift. 

This result i s cons istent with a recent compilation by 
iWatkins et~al"1 ()2009D who combined nine peculiar veloc- 
ity datasets measured using five different methods of dis- 
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Fig. 7. — Cosmological parameter likelihood surfaces in the flat toCDM model (one, two, and three sigma). The shaded contours display 
the K09 data set e (full sample) and the same contours are displayed in both panels for reference. The red cross indicates the point of 
maximum likelihood for the K09 data. In addition, the left panel shows an alternative analysis of the same data in which the CMB dipole 
correction is not applied (black curves). The right panel shows an alternative analysis where the both the CMB dipole and local dipole 
corrections are applied to the SN data. The point of maximum likelihood for each of these alternative analyses is indicated by the black 
diamonds. Only shifts perpendicular to the long axis of the contours are significant, because shifts along the long axis represent very small 
changes in \ 2 and are well constrained by other measurements (e.g. CMB and BAO). Correcting for the CMB dipole shifts the contours 
by about 15% of lcr, while the shift due to the local dipole is negligible. 



tance estimation (surface brightness fluctuations; fun- 
damental plane; type la supernovae; Tully-Fisher; and 
brightest cluster galaxies) and concluded that the bulk 
flow within a Gaussian window of radius 50/i _1 Mpc is 
407 ± 81km s- 1 toward I = 287° ± 9°, b = 8° ± 6°. They 
note that the magnitude of this flow is larger than would 
be predicted by standard cosmological models based on 
the best cosmolo gical parameter estimates from WMAP 
(jKomatsu et all 120091 ) . An even more significant de- 
viation from the predictions of ACDM was found by 
iKashlinskv et a l. (2008, 200jl), who measured the kine- 
matic Sunyaev-Zel'dovich effect of the CMB in the di- 
rection of known galaxy clusters. Their observations 
were at higher redshift (z ~ 0.1 or about 300/i _1 Mpc) 
and they found a considerably higher a mplitude bulk 
flow ( 600- 1 000km s" 1 ) than measured bv lWatkins et al.l 
(2009) and lHaugb0lle et al.l (|2007t ). but the dipole direc- 
tion was the same. As well as differing from the theoret- 
ical pred iction, this result does contradict the measure- 
ment by lHaugb0ll e et al.l (|2007) of a dipole amplitude 
that decreases with distance. 

Given the uncertainty in the redshift dependence of 
these local motions it is premature to apply any cor- 
rection before publishing observational data. Never- 
theless we want to estimate their effect on cosmolog- 
ical inferences. Since the low-redshift supernova data 
is the most influenced by redshift uncertainties, we es- 
timate local motions by a loc al dipole accord i ng to 
the consistent mea s ureme nts of iHaugbolle et aLl (2007) 
and iWatkins et"aT1 (|2009D at around z ~ 0.015. We 
then choose to u s e the redshift dependence measured by 
lHaugb0lle et al.l (|2007l ). in which the dipole magnitude 
decreases with distance. This relation is qualitatively 
what is predicted by standard ACDM models; however, 
quantitatively it does not drop as quickly as one ex- 
pects from simulations. Our choice of redshift depen- 
dence therefore lies between the theoretical predictions 
and the kinematic SZ results. The precise choice is not 



significant because it is the low-redshift points that have 
the largest dm/dz effect. 

The magnitude of the correction due to the local dipole 
is shown for the SDSS supernovae as a function of right 
ascension in Figure [5] and as a function of redshift in Fig- 
ure [SI The black points show the correction for the CMB 
dipole, ignoring the local dipole. The shaded points show 
the correction after the contribution from the local-dipole 
has been included to account for our lower velocity rel- 
ative to nearby galaxies than relative to the CMB. The 
lowest redshifts receive the largest local-dipole correc- 
tion, while the higher-redshifts tend towards the CMB 
correction. 

In Figure [7] we show the shift in cosmological parame- 
ter estimates that occurs when we ignore the CMB dipole 
correction and when we add the additional correction due 
to the local dipole. When we simulate the effect of these 
dipoles we change both the observed redshift and ob- 
served luminosity distance of the supernovae, according 
to Eq. 6 and 11 respectively. 

To make this comparison we have used the MLCS2k2 
version of data set 'e' from K09, which includes the 
new supernovae from the SDSS collaboration combined 
with the high-redshift ESSENCE, S NLS, and HST data 
along with the low-redshift sam ple (Ha muy et alJll996l : 
iRiess et al.l[l999t Uha et al.ll200l . 

The choice of dataset is important, so we also show in 
Figure[5]the effect of the dipole on data set 'b' from K09, 
which excludes the relatively isotropically distributed 
nearby sample and relies on the SDSS sample as the 
nearby anchor for the Hubble diagram. This makes it 
more sensitive to the dipole correction. 

It is clear from these figures that the CMB dipole cor- 
rection can be important. Neglecting it introduces a sys- 
tematic shift of Aw — 0.04 for data set 'b' when con- 
sidering the best fit w at a constant JIm ~ 0.3. This 
represents about 20% of one standard deviation at the 
current accuracy levels. This CMB correction is already 
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Fig. 8. — Similar to Fig. [3 but using dataset 'b' from K09, which includes the SDSS, ESSENCE, and SNLS data. This dataset is more 
sensitive to dipole corrections because it excludes the relatively isotropically distributed nearby sample and relies on the SDSS sample as 
the local anchor of the Hubble diagram. The best fit values are shown (as red crosses for K09, black diamonds for the two variations), but 
these are not particularly good indicators because the best fit in this case is on the edge of the parameter space explored. More indicative 
of the magnitude of the effect is the amount that the contours shift. The CMB dipole shifts the contours by about 0.2ct, corresponding to 
a Aio ~ 0.04 along the line of constant Q]^ = 0.3, while the shift due to the local dipole is small (Ami 5: 0.01). 



routinely applied to supernova data sets and it can be 
done to very high precision thanks to the accurate mea- 
surements of the CMB dipole. However, the characteris- 
tics of the local dipole are much more uncertain. Given 
the assumptions we have outlined here, Figure [7] shows 
the contribution from the local dipole is currently negligi- 
ble, giving a shift of Aw Z: 0.01. For the current sample, 
therefore, K09 are justified in correcting solely for the 
CMB dipole and ignoring local non-CMB contributions. 
This may not be true in the future when more data are 
available, particularly data at redshifts below 0.05. 

The lower sensitivity to dipole corrections shown in 
Figure [7] compared to Figure [5] demonstrates that choos- 
ing an isotropically distributed local supernova sample 
protects us, to a great extent, from systematic errors 
due to any unaccounted-for local dipole, because includ- 
ing all directions increases the scatter about the Hubble 
relation without adding bias. In the next section we will 
see that isotropic samples do not save us from the effects 
of higher-order motions. 

3.2. Impact of correlated velocities 

We use a model linear power spectrum and growth 
function based on a fiducial flat-ACDM cosmology with 
[h,Q m ,n b ,a 8 ,n] = [0.701,0.2792,0.046,0.817,0.96] to 
estimate the covariance matrix for the K09 sample, Cjj , 
as per Eq. [29] (see Fig. [4]). This encodes the likelihood 
that two supernovae will have correlated velocities based 
on their physical separation. 

Using this covariance matrix, rather than the usual 
uncorrelated error estimates, we re-fit our cosmological 
models. We fit a flat wCDM model, allowing the mat- 
ter density ft m and dark energy equation of state, w, 
to vary. (Note that we do not redo the velocity co- 
variance approximation for each different model, but the 
differences would be small.) To the supernova fits we 
add the same additional observations as K09, described 
in detail in their Section 8. Specifica lly, for Baryon 
Acou stic Oscillations (BAO) we use the ([Eiscnstci n et al.l 
120051) result that the derived distance parameter, A(z = 



0.35) = 0.469 ± 0-017 , and from the CMB we use the 
(|Komatsu et al.l I2009D result that the shift parameter 
R(z = 1100) = 1.710 ±0.019. 

The low-redshift cutoff is usually applied in order to re- 
move the effect of low-redshift peculiar velocities. Imple- 
menting a correlated velocity correction increases the er- 
ror bars on all the low-z supernovae relative to the high-z 
supernovae. This effectively down-weights the lower red- 
shift end of the Hubble diagram and thus has a similar 
effect to the low-redshift cutoff. 

We first consider only the low-redshift supernovae, for 
which this gives the largest effect, i.e. the Low-z+SDSS 
sample 'c' from K09. As done in K09, we apply a 
low-z cut of z cut = 0.02. In this case w increases by 
0.02 when the distance correlations are included. Using 
the larger sample 'd' from K09 that also includes the 
higher redshift ESSENCE and SDSS data, we find w in- 
creases by 0.014 due to the correlated velocities. These 
results should be compared to the uncertainty on w of 
±0.07(stat) ± O.ll(sys) reported by K09. Thus neglect- 
ing coherent velocities represents a potential systematic 
error on the best-fit value of w of up to about 2%, or 
about 13% of the current estimated systematic error 
budgetH When future supernova surveys achieve (su- 
pernova only) statistical error bars less than about 2%, 
this potential systematic error will need to be consid- 
ered carefully, especially for surveys with many nearby 
supernovae. Indeed the Carnegie supernova project have 
already found that the magnitude scatter in their sample 
of z < 0.08 supernovae is limited by peculiar velocities 
(see iFolatelli et aT]l20Tol Fig. 19). 

Figure [9] shows the effect of implementing a range of 
different low-redshift cutoffs on the supernova data, both 
for the original K09 MLCS2k2 data (solid lines) and for 
our version of that data with uncertainties corrected for 



28 Adding 0.014 to the systematic error budget represents an 
increase of 13% over the current 0.11 systematic uncertainty esti- 
mate. 
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Fig. 9. — Effect of low-z cut on SN data, where the solid lines 
show the best fit equation of state obtained while neglecting corre- 
lated errors and the dashed lines show the same when correlations 
have been taken into account statistically (model assumes flatness 
with w and Vl m the only free parameters, c.f. K09 Fig. 21). The 
upper and lower sets of lines show how much the results differ 
when you respectively include or exclude high-redshift data. Up- 
per (grey) curves are for the K09 data set 'd', with Nearby, SDSS, 
ESSFNCE, and SNLS supernova samples. Lower (black) curves 
are for the K09 data set 'c', with only the Nearby and SDSS data 
included. K09 use a low-redshift cutoff of 0.02 (crosses) with sys- 
tematic uncertainties (calculated for the flat-ACDM model with 
the MLCS light-curve fitter) of ±0.11 for set 'd' and +0.10 - 0.33 
for set 'c' (the greater systematic uncertainty in the lower direction 
for the set 'c' arises primarily due to uncertainties in the rest-frame 
U band). So the offset seen here between these two data set com- 
binations is within the uncertainties. The effect of correlations is 
currently smaller than the other systematic uncertainties consid- 
ered in K09, but will be important for attempts to measure w to 
better than 3%. Raising the low-redshift cutoff to 0.025 is sufficient 
to remove the expected effect of correlated supernova motions. 

correlated motion (dashed lines) q We plot the best fit 
value of w derived for a flat model with w and Sl m as 
free parameters. 

The correlation-corrected result (dashed line) can be 
matched by implementing a larger low-z cut on the 
uncorrected data (e.g. the correlation-corrected result 
with z cu t ~ 0.015 matches the uncorrected data with 
z cut pa 0.017). It is also evident that the effect of the low- 
z cut on the data with correlated errors is smaller than on 
the data set with uncorrelated errors. Both of these fea- 
tures are as expected, because some of the low-z cut was 
already effectively implemented by the down-weighting 
due to correlations. We also note that the SDSS su- 
pernovae are much less prone to correlations than the 
Nearby sample, simply due to their greater physical dis- 
tances between the supernovae. SDSS supernova there- 
fore provide a larger improvement to the low-redshift an- 
chor of the magnitude-redshift diagram than one might 
naively expect. 

This analysis demonstrates that if the velocity covari- 
ance matrix had been used in K09 then the need for a 
low-z cut would have been diminished. Using the covari- 
ance matrix for the peculiar velocity uncertainties should 
be an optimal statistical treatment of the supernova data. 
It automatically includes the effect of monopole uncer- 
tainties and dipole uncertainties as well as the higher- 
order correlated motions. Although slightly more com- 

29 Note that K09 also uses the SALT II light-curve fitter and 
the results differ. We do not debate the merits of different light- 
curve fitters here, our qualitative results are relevant whichever 
light-curve fitter is used. 



plicated to implement than a simple low-z cut, the results 
are more robust. 

Computing the full correlation matrix does have some 
disadvantages, primarily because it is model dependent. 
Our calculation of correlations has been made in a fidu- 
cial ACDM model, so it is not strictly self-consistent to 
use these correlations to test other models. However, 
this is mitigated by the fact that the majority of the co- 
variance signal comes from low redshifts, and most viable 
models for the universe have to agree fairly closely on the 
evolution and growth of structure in the local universe 
in order to match observations. To check the differences 
are negligible the correlations can be self-consistently re- 
derived for each cosmological being fitted. 

As future surveys with many more supernovae at- 
tempt to obtain percent-level accuracy on the value of 
the equation-of-state parameter the effect of correlations 
will become ever more important. At the very least, ne- 
glecting correlations under-estimates the uncertainty on 
our cosmological inferences, and in the worst-case sce- 
nario can bias the values of cosmological parameters we 
derive. 

4. IMPACT OF LOCAL UNDER-DENSITY 

Peculiar velocities are not the only source of system- 
atic redshift effects on nearby galaxies. The potential 
presence of a systematic redshift due to our possible 
position at the center of a local underdensity is of 
particular interest to supernova cosmology because the 
discovery of the acceleration of the universe is founded 
on the observation that high-redshift type la supernovae 
appear to be mor e distant than ex p ected in a decel- 
eratin g universe (jRiess et al.l 119981 : iPerlmutter et al.l 
1999). This observation has incited intense scrutiny 
of the potential for a local underdensity, known as 
a "H ubble bubble", tha t may be influencing our re- 
sults (jZehavi et al.l 119981: Uha et al.l 120061: iConlev et all 
l2007h . but the massive size needed to explain away 
the acceleratio n (lAlnes et al.ll2006t lEnavist fc Mattssonl 
1 20071: lEnavistl 120081: IGarcia-Bellido fc Haugbolld 120081: 
Garcfa-Bellido fe Haugbolld I2008L 120091: Izibin et al l 
20081: iClifton et all ^ 120081: I Alexander et al.l 120091: 
Marra fc Paakkonenl 120101 ) and the need to be 
almost at the cent er of suc h an underdensity 
(lAlnes fc Amarzguioul 120061: iCaldwell fc Stebbinsl 
120081 : iBlomqvist fe Mortselll l2010h makes such a model 
very contrived. 

Observations of local large scale stru cture indicate that 
we do sit in a local un derdensity (e.g. lGeller et al.|[l997l : 
IGottlober et al.l 120101) , albeit much smaller than that 
needed to explain a cosmological constant (on the or- 
der of 100 Mpc, as opposed to 1 Gpc). In Appendix C 
we use numerical simulations to show that underdensi- 
ties of this size are the size and depth of typical density 
fluctuations in a ACDM model. It is therefore impor- 
tant to investigate the impact that these realistic density 
fluctuations would have on our derivation of cosmological 
p arameters. 

iSinclair et al.l (|2010T ) used the Lemaitre-Tolman-Bondi 
(LTB) model to simulate a local underde nsity, as out- 
lined in IGarcia-Bellido fe Haugb0lld ()2008D and fit for a 
homogeneous cosmological model. They found that if we 
live in a 30% underdensity of scale 70/i _1 Mpc, then as- 
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Fig. 10. — This figure demonstrates the effect of correcting the SN data for a hypothetical 'Hubble Bubble', in this case a Gaussian 
undcrdcnsity of 30% on a scale of 70h~ x Mpc. In the left panel the shaded contours display the K09 data set e (with the maximum likelihood 
indicated by the red cross) and black outline indicates the results for the same data corrected for a local void (with the maximum likelihood 
indicated by a black diamond). In the left panel the low-redshift cutoff was the standard z cu t = 0.02. The right panel demonstrates the 
effect of increasing the low-redshift cutoff from z cu t = 0.02 to z cu t = 0.04. The analysis is identical to the left panel except that in both 
the homogeneous case (shaded contours) and the putative 'Hubble bubble' case (black outline), the low-redshift cutoff was z cu t = 0.04. 
Although increasing the low-redshift cut reduces the susceptibility of the data to local density fluctuations, dropping the low-redshift data 
changes the best fit cosmology by more than the local void would because of the weakened constraints. See also how the best fit w changes 
with changing low-z cut in Fig. [5] 



suming the universe is homogeneous could dupe us into 
believing that w is far more phantom- like (< —1), and 
the density of dark energy far more significant, than they 
really are (up to a 10% error in w in the worst-case sce- 
nario, even with a low-redshift cut that excises the local 
underdensity) . 

Given that these discrepancies remain significant it is 
worth investigating the potential impact on recent data 
sets and whether a higher low-z cutoff mi ght be worth- 
while. To that end we here apply the iSinclair et al.l 
(|2010D result to the SDSS data explicitly. At each 
observed redshift we alter the SDSS data for the ef- 
fects of a hypothetical void of the type described above 
(r = 70/i _t Mpc, with 5 = -0.3). The difference is 
shown in the left panel of Fig. [TUJ in which the standard 
low-z cut of z cut = 0.02 has been used. 

Although in Fig. HOl there appears to be a large shift in 
the best fit parameters, from (fi m , wq) = (0.40, —1.00) to 
(0.37, —0.86), the final cosmological parameters are not 
susceptible to this full discrepancy because most of the 
variation is along the long-axis of the contours, which is 
well constrained by other observations such as the CMB 
and BAO. The direction that the supernovae constrain 
most tightly is only changed a very small amount by the 
presence of a local void. 

It could nevertheless be argued that a higher low-z 
cutoff may be necessary to avoid the effects of local in- 
homogeneities. However, by excluding low redshift data 
one sacrifices constraining power. In the right panel of 
Fig. [10] we show the result of increasing the low-z cutoff 
to 

•^cut — 0.04. Although the effect of the void is now 
negligible (shaded and black contours overlap) the con- 
tours have shifted further due to the loss of constraining 
power than they did due to the hypothetical void. We 
therefore conclude that increasing the redshift cutoff is 
counter-productive with the size of the low-z SN sample 
used in K09. 



This analysis strengthens the argument for using the 
covariance matrix approach to down- weighting low-z su- 
pernovae since it inherently takes into account potential 
monopole velocities. 



5. CONCLUSIONS 

From this study we can conclude that the cosmolog- 
ical results derived by the SDSS SN survey are robust 
to peculiar velocity systematics. The local dipole repre- 
sents a negligible addition to the CMB dipole correction 
that has already been implemented. Future surveys with 
many nearby supernovae may need to take it into ac- 
count, but we note that an isotropically distributed local 
supernova sample would shield us, to a great extent, from 
systematic errors due to the local dipole. 

Neglecting correlated peculiar velocities can cause an 
error in the best-fit value of w, which in the current 
sample underestimates w by about 2%. It also causes 
us to overestimate the precision of our measurement. 
As future surveys aim for percent-level accuracy on the 
value of the equation-of-state parameter, the importance 
of correlations between the peculiar velocities of super- 
novae will increase. Here we treated them in a statistical 
sense, but it may be possible in the future to correct the 
supernova velocities for measured local flows. A future 
method of testing cosmological parameters will be to use 
the peculiar velocities as signal rather than noise, and 
generate a peculiar velocity power spectrum to compare 
against cosmological models. In the meantime, we find 
that accounting for peculiar velocities by using a covari- 
ance matrix for the correlated errors is a more robust 
way to down- weight low-redshift supernovae than apply- 
ing a sharp low-redshift cut. Doing so does not degrade 
the uncertainties in w, despite the down- weighting of the 
signal, because one can include more low-redshift super- 
novae in the overall fit. 
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Finally, we used n-body simulations to gauge the likely 
distribution of local under- and over-densities and found 
that a density fluctuation of 30% from the mean cos- 
mological density, out to a range of 70/i _1 Mpc, is rea- 
sonable given the expectations of concordance ACDM. 
A density fluctuation of this size can have a significant 
impact on the cosmological parameters we derive. The 
worst of those systematic errors can be avoided by down- 
weighting nearby supernovae and we demonstrated that 
the currently used z cut = 0.02 is well justified. However, 
we advocate including the velocity covariance directly 
in one's likelihood analysis as a more systematic way to 
down- weight the low-z SNe. 

In summary, peculiar motions and gravitational ef- 
fects due to inhomogeneities have systematic effects on 
the measurement of cosmological parameters using lu- 
minosity distance indicators such as type la supernovae. 
These effects will become significant for the next genera- 
tion of surveys and here we have suggested a covariance- 
matrix approach to correcting for them statistically, 
based on the expected correlation between the motion 
of the sources. This improves on the usual low-redshift 
cut approach, and we have provided code that generates 
the covariance-matrix for any sample of supernovae, to 
make this technique easy to implement. 
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APPENDIX 

APPENDIX A: TREATMENT OF RANDOM PECULIAR VELOCITY CONTRIBUTIONS 

The motion of distant supernovae and their host galaxies imprints a peculiar velocity error that is primarily random 
(as opposed to our own motion, on which Sect. l2,ll concentrates). That peculiar velocity dispersion crP cc ~ 300kms _1 
gives a redshift error of a pec = a pec /c (or the special relativistic formula if the peculiar velocity was higher). The 
measured redshift, z, is a combination of the recession and peculiar velocity contributions according to (1 + z) = 
(1 + z)(l + z pcc ), where z and z pec are the recession and peculiar velocity contributions to the redshift, respectively. 
Differentiating this expression to calculate the error contribution from peculiar velocities gives 

a z = (1 + z)<jP cc + (1 + z po >,-. (Al) 

We can take the error in recession velocity to be zero, a z = 0, so the uncertainty we need to add to our redshifts to 
account for peculiar velocities is a z = (1 + z )a pec . 

Previous analyses, including for example iDavis et al.l (|2007f) . used a z = a pcc and so slightly underestimated the 
contribution from peculiar velocities at high redshifts. Formally, the uncertainty at z = 1 should have been double 
what was used, but since the proportional contribution from peculiar velocities still decreases with redshift, this only 
corresponds to an error of 0.26%, as opposed to 0.13%, and the difference is negligible for cosmology. 

We convert a z into an approximate magnitude uncertainty, er^f , using the magnitude- redshift relation, and combine 
it in quadrature with the uncertainty in the measured magnitude, <7™ cas , and the intrinsic magnitude dispersion, <j™ , 
of the supernovae. 

30 http://healpix.jpl.nasa.gov 
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FIG. 11 Examples of the conver- 
sion from redshift uncertainty 
to magnitude uncertainty. A 
peculiar velocity uncertainty of 
crfj cc =300kms~ 1 corresponds 
to o-p cc =0.001. This value con- 
verts to a large magnitude un- 
certainty at low redshift, where 
the slope of the magnitude- 
redshift diagram is steep, but a 
smaller magnitude uncertainty 
at high redshift. Different fidu- 
cial models give slightly dif- 
ferent conversions between red- 
shift and magnitude uncertain- 
ties, but the difference is negligi- 
ble for cosmological inferences. 
Here the emp ty model conver- 
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The distance modulus has the form 



/i = 51og 10 (d L ) 



In [x(l + z)\ , 



(A2) 



where x = RoX = 
magnitude error of 



ln(10) 

: J* H(z)^ 1 dz is the comoving distance. Therefore an error in the redshift corresponds to a 

5 I 



ln(10) 



1 



xH(z) 



(A3) 



Although a fiducial cosmology is used for this calculation - often taken to be ACDM with Om ~ 0.3 and J7a ~ 0.7 - 
differences from the derived cosmology are small and have negligible impact on cosmology fits. Note that K09 use the 
empty universe as their fiducial cosmology for error calculations, and approximate the empty universe case, in which 
H(z) = H (l + z) and x = cln(l + z)/H , by 



ln(10) 



1 + z 



z(l + 1/2) 



(A4) 



The difference between this approximation and Eq. IA3I is shown in Fig. 111! 

In the non-flat case x should be replaced with RoSk(x) m t ne equation for fx, or RqT^x) m the equation for er^, 
where Sk — sin or sinh in the closed and open cases, respectively, and Tk — tan or tanh. 

APPENDIX B: THE EQUIVALENCE OF THE OBSERVER-CENTRIC AND SEPARATION-CENTRIC EXPRESSIONS FOR 

THE MAGNITUDE COVARIANCE MATRIX 

Here we derive the expressions for Clf in Eq. [29] and [30] These derivations are valid in the flat universe case. It 
comes down to evaluating the two point velocity correlation ^| = (( v i ' xi)( v 2 ■ £2)), where 1 and 2 label the two 
SNe in question, since 
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This is actually an old subject (see e.g. iGorsk i 1988). One reason we go over the derivation here is that errors have 
crept into some recent literature, as pointed out by iGordon et~aT1 (|2007T) . It is also useful to see how two completely 
different looki ng: expressions, i.e. Eq. 1291 and 1301 ar e actually equivalent. Err ors have occurred in some recent versions 
of Eq. [29] (e.g. iHernandez-Monteagudo et alJl2005 iCoorav fc Caldwelll[2006l ). 
Using linear theory, it can be shown that 

j^k- 2 P(k) z=Q (k • Xl )(k • x 2 )e- k -( Xl - X2 > (B2) 

where Xi and x 2 are the comoving positions of the two SNe in question, Xi and x 2 are the unit vectors pointing in 
these directions, P(k) z= o is the mass power spectrum today, and D[ and D' 2 are the derivatives of the growth factor 
with respect to conformal time at the two redshifts of interest. 
An observer-centric approach is to use 

k • x 2 e* k ^ = 47r^- 1 j^;b)^/ m (k)W*2) (B3) 

where ji is the spherical Bessel function, j' t is its derivative (with respect to its argument, not conformal time), and 
Yarn's are the spherical harmonics. Performing the integral over k in Eq. lB2[ and using J dflfcYf (\s.)Yii m i (k) = Su'S m m' 
and Ve(xi ■ x 2 ) = 4n/(2£ + 1) J2 m ^m( Xl )^m( x 2) , it is straightforward to show that 

/dk 
-^P(fc) z=0 J2( 2i + 1)J/(*Xx)jK*X2)^(xi • x 2 ) (B4) 

from which Eq. [29] can be obtained (see lHui fc Gr eene 20061 for details). The above expression is observer-centric in 
the sense that one can easily read off from it the angular velocity power spectrum as seen by the observer, 

—P(k) z ^MkXi)A(kX2). (B5) 

Here, 1 and 2 can refer to the same re dshif t, or two different redshifts. 

A different approach to reducing Eq. IB2I is to first note that by symmetry arguments (Gorski 1988), 

MxOt^xa)) = P(r) - ^(r)]fifj + E(r)*«, (B6) 

where i and j here, unlike in the rest of the paper, label the spatial directions rather than the SNe, r is the comving 
separation between points 1 and 2, and f is the associated unit vector. Suppose f points in the z direction, then the 
above matrix is diagonal, with diagonal entries S,S,II i.e. £ is the perpendicular velocity correlation and II is the 
parallel velocity correlation. Here, parallel and perpendicular are defined by the separation vector between the two 
SNe (hence a separation-centric approach). From this matrix, one can deduce that 

^2 Cl - (xi ■ f)(x 2 • f)H(r) + • x 2 - (* x ■ f)(x 2 ■ f)]S(r) (B7) 

where [xi • x 2 — (x i ■ f ) (x 2 • f )] can be written as sm8i sin6*2 if xi • f = cos^i and x 2 • f = cos6*2 . Comparing this 
expression with Eq. IB2| one can see that, 

/r1 3 k 
^fc- 2 P(fc) 2=0 (k-f)V k - r . (B8) 



Using 

e ik r _ 



e* r = y (21 + l)i e Je (kr)Vi(k ■ f) (B9) 



and integrating over k (choosing f to lie in the z direction for instance), one can see that only I = 2 and I = survives. 
Finally, using j 2 — 3j%/x — jo, one obtains, 



H(r) = D[D' 2 J ^P{k) z =o 



Jo(kr) - 



kr 



(BIO) 



The perpendicular counterpart can be similarly obtained from, 



d 3 k I'ik). , - , ;k ,. 



S(r) = D[D' 2 j {27t)3 " k2 (k ■ *ye"', (Bll) 



with x pointing in the x direction while f points in the z direction. A few manipulations yield, 

m=D' l D' 2 J^P{k)^^-, (B12) 
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Fig. 12. — Left: The average monopole (blue) and dipole (green) in the velocity field as extracted from a large N-body simulation with 
best-fit WMAP5 cosmological parameters. The shaded areas indicate the cosmic variance. The dashed line is the best-fit model, Eq. JUT}. 
Right: To what distance can the direction of the local dipole be extrapolated? This plot shows the relative direction of the dipole as 
measured by observers at two reference redshifts zo = 0.015, 0.2 (blue and green respectively). When the change reaches 90° (dashed lines), 
which occurs at 290=0. 043 and zgo={0. 132, 0.294}, the dipole to a shell at that radius bears no correlation with the dipole in shells very 
close to the observer. The shaded area is the cosmic variance. As expected, the size of the region in which the directions of the dipole are 
aligned increases with redshift, from AZ90 ~0.03 to Azgo ~0.08 for the two cases shown here. 



reproducing the results of IGorsk i (1988) and giving our Eg. I3U1 

APPENDIX C: MODELLING THE SIZE OF DENSITY FLUCTUATIONS 

Once we have derived a model (such as ACDM) from the observational data we can ask whether the density 
fluctuations predicted in this model are consistent with those observed. We can also address whether our treatment of 
the dipole is justified and whether it is likely that an undiagnosed monopole term could be biasing our results. This 
checks the internal consistency of the model as well as testing for biases our assumptions may impose on our results. 

To this end we performed a large scale d ark matter N-body simulation (Lb O x=2048 Mpc h~ l , A r part =1024 3 , z s tart=49) 
using the Gadget2 code ( |Springell |200 5j) with best-fit WM AP5 cosmological parameters {Q m ,Q\, h,w,n s , a s } = 
{0.2792, 0.7208, 0.701, -1, 0.96, 0.817}. (|Komatsu et al.ll2009f ). Using 2000 observers placed at random, but weighted 
by mass, we calculated the average magnitude of the monopole and the dipole, together with the cosmic variance 
of each (see Fig. [T2")) . We define the magnitude of the monopole to be a(M) = yf^{\M\), where a(M) is the root- 
mean-square of the signed monopole, and the signed monopole is simply the mean velocity (either towards or away) 
of matter in a shell of a particular distance. 

While the analysis is done at redshift zero using a single data snapshot, the velocities have been corrected using 
linear theory, so the velocities in Fig. Q2] are in the lightcone. The correction from translating the velocities to the 
lightcone is minor, at maximum 1%. 

In order to calculate the mean size of the dipole and how it varies with distance (redshift) we sliced the resulting 
simulation into shells with a range of radii between 10 h~ 1 Mpc < r < 1000 /i -1 Mpc around each of the random 
observers. We then measured the mean motion of these shells to calculate the monopole and dipole the central 
observer would see for sources at that distance. 

We have done the analysis both with shells with a thickness of 10 ft, -1 Mpc and of 1 /i -1 Mpc, and confirmed that 
the results do not depend on the shell thickness, except at the very lowest redshifts where the thickness becomes 
comparable to the radius of the shell. The magnitude of the dipole is reasonably well described by the simple model 

v d = (507 ± 51) - (65 ± 8) In ( ) kms" 1 , (CI) 

= (-8 ± 12) - (63 ± 7) ln(z) km s" 1 , 

where the differences in the fits using either comoving distance (i?ox) or redshifts are due to the slightly nonlinear 
conversion between the two at larger distances. 

The mean magnitude of the dipole at low redshifts (z ~ 0.01) is approximately 300 ± 100 kms -1 . This is consistent 
with the average random peculiar velocity uncertainty we assume for supernovae. 

The mean absolute magnitude of the monopole, which is also plotted in Fig. I12[ is smaller than the dipole but still 
significant, on the order of 100 kms -1 . We investigate in Section 0] the impact this mean monopole would have on our 
cosmological inferences. 

The direction of the dipole of the local velocity field is only known at very low redshifts. To test how well this 
knowledge can be extrapolated to higher redshifts, Fig. Q2] shows how the direction of the dipole in the simulation 
changes as a function of redshift. The direction is measured with respect to the direction of the dipole at two reference 
redshifts z = 0.015,0.2. These redshifts correspond respectively to the redshift of the currently available local dipole 
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measurement, and a redshift representative of the SDSS supernovae. 

A naive expectation would be that the average local dipole should decrease with redshift until we reach sources that 
are too distant to share any significant common source of gravitational attraction with us. At that point the sources 
should be on-average at rest with respect to the CMB, and therefore our dipole direction with respect to those sources, 
if we do not correct for our local velocity, should be simply the direction of the CMB dipole. When we do correct for 
the local velocity (as done in Figure H"2j) . while the amplitude of the local dipole decreases, the direction still changes 
at higher redshifts, and it only makes sense to extrapolate the currently known dipole direction out to z « 0.045. 

This result is interesting seen in light of observationa l results pointing towards a coherent dipole direction out to at 
least 300 /i~ 1 Mpc (i.e. z ~ 0.1, Kasminsky_et_alJ[2008) , since not only the magnitude of the observed dipole velocity 
(1600±500kms _1 at z ~ 0.03 and 850 ± 250kms~ 1 at z ~ 0.1) but also the constancy of the the direction of the 
dipole is surprising (~ 3cr deviation) when interpreted in the framework of the ACDM cosmology. 
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